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1.  Introduction 

The  prediction  capability  of  a turbulence  model  depends  on  how  effec- 
tively one  can  prescribe  the  Reynolds  stress  distribution  in  closing  the  system 
of  equations.  The  simplest  and  most  widely  used  has  been  the  Boussinesq 
treatment  of  the  Reynolds  stresses.  As  is  now  well  known,  Boussinesq  hypo- 
thesis holds  only  when  the  strain  rates  are  fairly  small.  The  main  reason 
being  that  by  construction  the  Boussinesq  formula  implies  that  the  principal 
axes  of  the  Reynolds  stress  tensor  are  parallel  to  the  principal  axes  of  the 
strain  rate  tensor  so  that  any  change  in  the  strain  rate  is  directly  felt  in 
the  stresses.  This  instantaneous  change  of  the  Reynolds  stresses  with 
the  strain  rates  is  not  supported  by  the  experimental  observations,  because 
the  Reynolds  stresses  being  due  to  the  vorticity  fluctuations  require  some 
time  to  adjust  to  the  new  strain  rates.  To  overcome  these  short  comings, 
one  must  cither  abandon  the  Boussinesq  hypothesis  altogether  and  solve  the  six 
Reynolds  transport  equations  which  is  costly  in  terms  of  computer  time  or 
improve  upon  the  hypothesis  itself. 

In  this  paper  we  follow  a recent  analysis  of  Rodi^  to  construct  an  im- 
proved second-order  version  of  the  Boussinesq  hypothesis.  An  algebraic  rela- 
tion for  the  turbulent  stresses  has  been  obtained  through  a consideration  of 
the  transport  equations  of  the  Reynolds  stresses.  Consequently,  the  resulting 
relation  has  tlie  necessary  influence  of  the  convective  and  diffusive  trans- 
port effects  of  a turbulence  stress  field. 


This  paper  is  an  outgrowtii  of  a current  research  supported  by  the  United 
.Slates  Air  Force  Office  of  Scientific  Researcli,  Grant  No.  AFOSR-7G-2922 . 
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In  this  paper  an  algebraic  relation  for  the  Reynolds  Stresses  has  been  obtained 
through  a consideration  of  the  transport  equations  of  the  Reynolds  Stresses. 
ITnis  analysis  provides  a second-order  approxicat ion  to  the  Boussinesq  eddy 
viscosity  hypothesis.  The  basic  assumption  of  the  analysis  is  that  the 
derivatives  of  the  ratio  of  the  Reynolds  Stresses  and  energy  is  small  in 
comparison  with  the  other  terms. 
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2.  Analysis 


The  transport  equations  of  the  Reynolds  stresses  (-u^u^)and  the  equation 
of  turbulence  energy  (e  = ~ incompressible  flow  respectively  are 


^ = P + D - e 
at 


(1) 

(2) 


where  — is  the  substantive  derivative  based  on  the  mean  velocity  components 
dt 

^ij*  ■ ’ ^i  ‘ production,  diffusion  and  dissipation 

of  the  Reynolds  stresses,  is  the  pressure-strain  correlation,  while  P, 

D and  c respectively  are  the  production,  diffusion  and  dissipation  of  the 

turbulence  energy.  In  this  paper  we  have  utilized  the  modeling  of  the  terns 

0 . D , E . and  D as  reported  in  references  1 and  2,  which  on  using  the 

'<ij’  ij’  i3 

summation  convention  on  repeated  indices  are 

c,  e 


- Y <3  ^ hj  - ’13>  hj  - 

, 3t . . c e 9t . . 

D . . = r — ( ^ T,  ^ 

ij  9Xj^  3x^  E k<-  ?x. 


(3) 

(4) 


‘'ij  3 ^"ij 

n = _i_  (^JL  + ^ JL) 

3Xj^  9xj^  E ^kC  9x. 


(5) 

(6) 


vlii're  y and  c^  are  empirical  constants,  and  v the  kinematic  viscosity 


'l!ie  ti'im.s  P..,  p .and  c are 


ij 


?U, 


?u. 


P = -(t  — 3-  + 1") 

ij  ''  ik  9x^  jk  9x,  ^ , 

j 9U 

' .U  , 
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h',  ■ 
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5V  , , 


r’y 


Introducing  the  notation 


3.1 


’ij/ 


(1 


7 


(7) 

(b) 

(V) 

(10) 


and  arranging  terms  in  D , we  have 

i . - DT, , . 

D. , * v-^  (e  + v-~  + DT.  . 

Ij  ij 


c^(e)2 


T _Jj. 


c e - DT. . 

c kl  D,^  Dxj^ 


Introducing  (11),  the  identity 

dT_  dT..  ,- 

_li  = -_iL  + T ^ 
dt  dt  Ij  dt 

and  Eq.  (2)  in  Eq.  (1)  and  neglecting  the  derivatives  of  in  comparison 
with  the  other  terms,  we  obtain 

Ty(P-c)  - Pj.  Qjj  - Cjj  (12) 

On  substituting  (3)  and  (5)  in  (12)  we  obtain 

^ij  “ I " I 


where 


Y = 1 - Y and  d,  = c,  -1 
o 11 


We  now  introduce  the  following  notation 

, , 9U. 

■ 2"l“l  • “ij  ■ 6 ^ <“> 

J 

where  9 is  the  vorticity-density  and  is  the  fluctuating  vorticity  compo- 
nent. It  follows  directly  from  the  Kolmogorov-Saf fman  equation  of  energy 
(Ref.  3)  that  the  dissipation  of  energy  c is  given  by 

c = e6  (15) 

I'siiig  (14)  and  (15)  in  (7)  and  (8)  we  get 

■ -(’ik  “jk  ^ hk  “ik> 

V/c  = -Tj..  Mj^  . (17) 

If  wc  now  substitute  (16)  and  (17)  in  (13)  then  we  get  a system  of 
iioiil  invar  simultaneous  algebraic  equations  for  the  determination  of 
R.hH^  in  his  derivation  did  not  use  the  expansion  (17),  but  retained  P/t 


as  a parameter  and  solved  (13)  for  j ■ Siiieo  P/c  contains  all  T^^'s, 

We  lollow  an  approach  different  from  Roili,  which  in  the  first  place  cstahli- 


sho<;  the  validity  of  the  Boussinesq  hypothesis,  and  In  the  second  place 
yields  an  improved  algebraic  relation  for  j • 

Equation  (13)  can  be  written  in  various  iterative  forms,  however,  the 
following  form  is  chosen  because  it  yields  various  approximations  in  a direct 


fashion. 


. _ (n+1)  (n)  (n) 

‘^l  \ji 


ij  W 


' I "jk  - "jk"'  “ik  ^ I "J"’  "k;,> 

where  n is  the  iteration  index. 

For  the  zeroeth  approximation  we  take 

T (0)  = 2 , 

'ij  3 ^ij 

and  by  using  the  continuity  equation  = 0,  we  get 

"if’ 

where 

2 

a = “ constant.  (20) 

o jdj^ 

Equation  (19)  is  the  usual  Boussinesq  formulation  which  has  also  been 

4 3 

used  among  others  by  Kolmogorov  and  Saffman  . It  is  only  a first  approxi- 
mation and  is  expected  to  hold  in  situations  where  the  mean  flow  is  not  chang- 
ing rapidly. 

To  obtain  the  second  approixmation,  we  introduce  (19)  in  (18)  and  neglect 


terms  of  the  third-order  in  M. . to  have 

ij 

(0)2 

* |,.‘r<Mj,  + + (M.,  + - I + M 1 (21) 

■ c ili.'ii  (21)  provides  the  second  approximation  to  the  Reynolds  stresses  ar..i 
• ‘."X'Vv- tod  to  hold  from  low  to  moderate  variations  of  the  strain  rates. 

By  lollowing  a pliilosopliically  different  approach  Saffman'  has  also 


.in  cxpressicMi  similar  to  (21)  wliich  in  our  notation  Is 


(22) 


’ '1  ‘u  ■ "jt’ 

-|'V^''kj>\i*<"u*V\j)  > 

where  u and  X are  constants.  Comparing  (21)  and  (22)  we  find  that  though  both 
are  in  dimensional  agreement,  they  diffei/  in  their  last  terms.  It  must  be  noted 

that  (21)  is  a consequence  of  the  complete  Navier-Stokes  equations  while 
(22)  is  an  attempt  at  finding  a relaxation  model  to  overcome  the  diffici- 
encies  of  the  first  approximation. 

A comparison  of  (21)  and  (22)  yields  the  values  of  the  constant  appearing 

in  (21).  Thus  we  have  i 

a = a = 0. 3 
o 

3 

The  value  a^=  0.3  has  been  used  by  Saffnan  and  also  by  Pope  and  Khite- 

law^,  but  from  (20)  we  find  that  the  value  =0.3  is  not  consistent  with  the 

o 

values  7^  = 0.4  and  dj^  = 0.5  as  proposed  in  Ref.  2.  However,  if  we  take 

2 6 

the  value  d,  = 1.86  - - 1)  as  mentioned  by  Rotta  and  =0.3  then 

i.  . / o 

~ 0.33,  which  is  near  to  the  value  used  by  Saffman.  For  the  sake  of 

definiteness  we  therefore  select  the  value  a •=  0.3  in  Eq . (21). 

o 

lor  two-dimensional  mean  flow  Eq . (21)  yields  the  following  expressions 
lor  rill'  normal  and  shear  stresses: 


’ll 

(2) 

_ T (1) 

11 

o 

o 

[ 2:-r 
‘ 11 

+ 

3(Mj2  + ■ 

■ ^‘\2  -'2i3^^ 

(23a) 

(2) 

.r  (1) 
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■> 

3(>!„  + - 

^21 

c T 

+ a 

o 

‘ 11 

+ 

■ ^‘‘12  •':i^^^ 

(23b) 

(2) 

T (0) 

4 

o 

’.1? 

33 

- a 
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(23c) 
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(2) 

_ T 
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ir  i; 
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(23d) 
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On  the  otlier  hand  Saffnan's  equations  (22)  yields 

^11  “ ^11  ^^^12  ”21^  ^^12  ~ ^21^ 


>.(M^2  + -^21^  ^'’12  “ ^21^ 


, (2) 

= T 

11 

^11 

, (2) 

= T 

22 

22 

, (2) 

_ T (0) 

33 

33 

, (2) 

_ T 

12 

- ^12 

2>.Nii(Mi2  - M^,) 


(24a) 

(24b) 

(24c) 

(24d) 


where  A is  an  emperical  constant.  Thus  though  the  shear  stresses,  Eqs.  (23d) 
and  (24d),havc  the  same  distributions,  the  normal  stresses,  Eqs.  (23a-c)  and 
Eqs.  (24a-c),  have  entirely  different  distributions. 

In  the  wall  region  = >1^2  = ^21~  ^ 

“12  " i 

a 

o 


so  that  Eqs.  (23)  become 

(2)  2 2 

T , - — = 2a 

3 o 


T (2)  _ 2 _ 2 

22  3 o 

a,  (2)  2 2 

^33  3 " 

^12  "“o 


Equations  (24)  become 

2 = ,/,  2 
11  3 o 

22  ”3 

^ (2)  _ 2 


(2) 

T ^ ' ^ 

12  o 


(25a) 

(25b) 

(25c) 

(25d) 

(26a) 
(26  b) 

(26c) 

(26d) 


Safin;.’.:',  now  tabes  \ =>  .02  to  n..ilcii  tiio  ihrco  r.orr..il  stressc.s  with  the  i Mi'ci  I- 
nental  u.Ua  which  roughly  aio  in  the  r.itio  4:2:3.  N'umoiical  v.ilues  ba.sed  oii 
(25)  and  (26)  and  the  expcri:;icr.ial  valuo.s  as  quoted  in  P.ef.  7 are  tabulated 


below. 


Table  1 Comparison  of  the  Ncar-rWall  Data 


1 


I 


i 

I 


I 


L 


Present(a=0. 3) 

3 

Present  (a=0.  34)  Saffnian 

Reference  7 

, (2) 
11 

2 

3 

0.18 

0.23 

0.22 

0.32 

, (2) 
22 

2 

~ 3 

-0.09 

-0.12 

-0.22 

-0.18 

. (2) 
33 

2 

3 

-0.09 

-0.12 

0.0 

-0.10 

, (2) 
12 

-0.30 

-0.34 

-0.30 

-0.34 

A comparison  of  values  in  Table  1 shows  that  a^=0.34  in  the  present  model 
may  be  more  suitable  than  a^=0.3.  However,  any  adjustment  of  the  constant 
or  the  actual  prediction  capability  of  the  proposed  second-order  algebraic 
relation  (21)  can  be  ascertained  only  after  it  has  been  used  in  the  calcula- 
tion of  various  turbulent  flows. 

2 

Based  on  the  works  of  Launder  et  al  and  on  the  most  recent  review  by 
0 

Reynolds  it  is  possible  to  establish  in  advance  the  limitations  of  the  second 

approximation,  viz.  Eq.  (21).  For  example,  it  may  be  observed  that  in  the  wall 

(2)  (2) 

region  the  normal  stresses  and  are  equal  (Eqs.  25b, c).  This  result 

is  in  exact  conformity  with  Eqs.  (14)  of  Ref.  2 in  which  Launder  et  al  have 
shown  that  in  any  simpler  pressure-strain  hypothesis  there  is  no  direct 
production  of  and  T^^  and  therefore  they  tend  to  be  equal.  Thus  it  can 
roughly  be  stated  that  the  amount  of  approximation  involved  in  (21)  is  the  same 
as  obtaining  the  Reynolds  stresses  through  the  six  Reynolds  stress  transport 
equations  by  using  a simpler  pressure-strain  relation. 
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3.  Conclusion 


The  purpose  of  this  paper  has  been  to  denonstrate  that  the  Boussinesq 
eddy  viscosity  hypothesis  and  its  higher  approximations  are  a direct  conse- 
quence of  the  Reynolds  stress  transport  equations.  The  basic  assumption  of 

the  analysis  is  that  the  derivatives  of  the  terms  T.  . = '^ij  are  small  in 

^ e 

comparison  with  the  other  terms  in  the  rate  equation  for  • Because  of 
the  implicit  effects  of  the  convection  and  diffusion  in  the  second  approxi- 
mation, the  simple  algebraic  relation  (21)  is  expected  to  provide  a basis 
for  the  prediction  of  com.plcx  flows. 

It  is  important  to  mention  that  the  validity  of  Eq.  (13),  which  rests 
on  the  assumption  that  the  rate  of  change  of  be  small  in  comparison  with 
the  other  tencs  is  hot  new  either  in  the  paper  by  Rodi‘  or  in  the  present  one. 
This  idea  5 has  earlier  been  used  by  Donaldson'  who  called  it  as  the  "super- 
equilibrium."  limit.  As  an  application,  Donaldson  obtained  all  the  Reynclds 
stress  terms  algebraically  for  a line  vortex. 


The  authors  arc  grateful  to  t'ue  rov:vwer  for  bringing  this  to  our 


attention. 
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